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o 

O 1 We numerically investigate the Josephson transport through ferromagnetic insulators 

(FIs) by taking into account its band structure. By use of the recursive Green's function 
method, we found the formation of the tt junction in the case of the fully spin-polarized 
FI (FPFI), e.g., La2BaCu05. Moreover, the O-tt transition is induced by increasing the 
thickness of FPFI. On the other hand, Josephson current through the Eu chalcogenides 
shows the tt junction behavior in the case of the strong d-f hybridization between the 
conduction d and the localized / electrons of Eu. Such Fl-based Josephson junctions 
may become a element in the architecture of future quantum information devices. 



iy~v , Keywords: Josephson effects; tt junction; spin filter; spintronics; Green's function method; 

^"^ 1 quantum computer. 

o 

1. Introduction 

(^j ■ Ferromagnet-superconductor hybrid struct ures exhibit novel phenomena which have 

'V^ . been studied extensively in the recent yearsJJ^ These systems provide the possibility 

C^ I for a controlled study of coexistence and competition of the ferromagnetism and the 

superconductivity. One of the most interesting effects is the possibility of forming 
the TT Josephson junction in s uperconductor /ferromagnetic-metal/superconductor 
(S-FM-S) heterostructures.'^^ In a tt junction the ground-state phase difference 
between two coupled superconductors is tt instead of as in the ordinary junctions. 
The existence of the tt junction in S-FM-S systems has been confirmed in experiment 
by Ryanzanov et al.l^ and by Kontos et al.'-^ In terms of the Josephson relationship 
I J — Ic sin (j), where (j) is the phase difference between the two superconductor layers, 
a transition from the to tt states implies a change in sign of Ic from positive to 
negative. Physically, such a change in sign of Ic is a consequence of a phase change 
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in the pairing wave- function induced in the FM layer due to the proximity effect. 

Recently, quiet qubits consisting of a super conducting loop with a S-FM-S n 
junction have been theoretically proposed ! II Ifn the quiet qubits, a quantum two 
level system (qubits) is spontaneously generated and therefore it is expected to be 
robust to the decoherence by the fluctuation of the external magnetic field. From the 
viewpoint of the quantum dissipation, however, the structure of S-FM-S junctions 
is inherently identical with S-N-S junctions (N is a normal nonmagnetic metal). 
Therefore a gapless quasiparticle excita tion i n the FM layer is inevitable. This 
feature gives a strong Ohmic dissipatio n^ ^ I and the coherence time of S-FM-S 
quiet qubits is bound to be very short. 

On the other hand, as was predicted by Tanaka and Kashiwaya,^'^ the n junc- 
tion can be formed in Josephson jun ctions with ferr omagnetic insulators (FIs). By 
using the functional integral method! ^ II I I we have theoretically proposed 
a superconductor phase ' and flux type qubit a° l ^ l ^ l based on S-FI-S n junctions. 
Moreover we have showed that the effect of the dissipation due to t he q uasi-particle 
excitations on macroscopic quantum tunneling is negligibly smalL^!^ However, in 
above studies, we have used a very simple 5-function model as the FI barrier. There- 
fore, the correspondence between this toy mode and the actual band structure of 
FIs is unclear. In this paper, we will formulate a numerical calculation method for 
the Josephson current through FIs by taking into account the band structure of 
FIs. Then we will discuss the possibility of the formation of the 7r-coupling for the 
Josephson junction the two types of the FLi.e., the fully polarized FI (FPFI) and 
the Eu chalcogenides (e.g., EuO and EuS).!^ 

2. Energy band structure of ferromagnetic insulators 

The typical density of states of FPFI for each spin direction is show n in F ig. 1(a). 
One of the representative material of FPFI is undoped La2BaCu05! I I The ex- 
change splitting Vex is estimated to be 0.34 eV by the first-principle band calculation 
using the spin-polarized local density approximation.'^^ Since the exchange splitting 
is larger and the bands were originally half-filled, the system becomes an FI. 

On the other hand, recently spin filtering effect are intensively studied by use 
of the Eu chalcogenides. The schematic energy-band structure of the Eu 

chalcogenides is shown in Fig. 1(b). The Eu chalcogenides stand out among the 
FIs as ideal Heisenberg ferromagnets, with a high magnetic moment and a large 
exchange splitting of the conduction band for Eu 5d-electrons. Utilizing the exchange 
splitting (V^x) to filter spins, these materials produce a near-fuUy spin-polarized 
current when used as a tunnel barrier. Of the Eu chalcogenides, EuO has the largest 
V^y. and the highest Curie temperature (Tcurie ~ 69 K for bulk). 

In EuO, the large saturation magnetic moment 7/iB per Eu^"*" originates from 
the seven unpaired electrons localized at the 4/ band below the Fermi energy. 
Ferromagnetic order of the 4/ spins causes exchange splitting of the conduction 
5rf band, lowering (raising) the spin-up (-down) band symmetrically by V^^/2. A 
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(b) 




Fig. 1. The density of states for each spin direction for (a) the fully polarized ferromagnetic 
insulator and (b) the Eu-chalcogenides. 

large exchange splitting of 0.54 eV was determined by measuring the redshift of the 

l27l 
absorption edge in single crystals of EuO cooled below T'curie-'^^ 

When an ultrathin film of the Eu chalcogenides is used as the tunnel barrier 

between two metallic electrodes, the exchange splitting of the conduction band 

gives rise to a lower barrier height for spin-up electrons and a higher barrier height 

for spin-down elec trons. Because of the tunnel current depends exponentially on 

the barrier height j ^ I the tunneling probability for spin-up electrons is much 

greater than for spin-down electrons, leading to a highly spin-polarized current. 

This phenomenon is called the spin-filter effect. 



3. Numerical calculation of Josephson current 

In this section, we develop a numerical calculation method for the Josephson current 
of S-FI-S junctions. Let us consider the two-dimensional tight-binding model for a 
S-FI-S junction as shown in Fig. 2. The vector r ~ jx + my points to a lattice 
site, where x and y are unit vectors in the x and y directions, respectively. In the y 
direction, we apply the periodic boundary condition for the number of lattice sites 
being W. 

Electronic states in superconductor are described by the mean-field Hamiltonian 

"Hbcs =2^1^'' ^'''''' ^'■' ~ ^■^'^ 
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Fig. 2. A schematic figure of a Josephson junction through the ferromagnetic insulator on the 
tight- binding lattice. 

with Cr = (cr,t, Cr,^)^ where c\. ^ (c^ g.) is the creation (annihilation) operator of an 
electron at r with spin a = ( t or J, ), c means the transpose of c, and (Tq is 2 x 2 unit 
matrix. The Fermi energy /i is set to be 4i for superconductors. In superconductors, 
the hopping integral t is considered among nearest neighbor sites and we choose 
A = iA(T2, where A is the amplitude of the pair potential in the s-wave symmetry 
channel, and a^ is a Pauli matrix. 

We will consider two types of FI as a barrier of the Josephson junction. In the 
case of FPFI, the Hamiltonian is given by a single-band tight-binding model as 

^FPFI = -^ X! ^r.<y^r'.a " ^^{^i " M)c|.,t"^^,T ~ 51 ('^* ~ MFPFI + Vc^)c\j^Cr^^, 
r,r',(T r r 

(3) 

where Vex is the exchange splitting [see Fig. 1(a)]. If V^x > 4i, this Hamiltonian 
describes FPFI. The Fermi energy /ufpfi is set to be Vox/2 — 4t. 

On the other hand, in the case of the Eu chalcogenides, we have used a following 
d-f hamiltonian, 

HEc = Hd + nf + ndf, (4) 

Hd = -td ^ dl,^dr',a - ^{4:td - f^d)dl^^dr,^ - ^{4:td - fJ-d + Vcx)di^idr,i, 

r.r',(T r r 

(5) 



Udf = Vdf ^ (4,<T/r-,a + fr,adr.a) , 

where dj. ^ {fr.a) i^ ^^^ creation operator, td (tf) is the hopping integral and V^ 



(6) 
(7) 
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{VjjJ is the exchange splitting of d{f) electrons. The Fermi energy of d and / elec- 
trons is respectively given by ^d — —Qd and /i/ = 8t/ + gj , where gd (g/) is the 
energy gap of the d{f) band [see Fig. 1(b)]. The third term Hdf of the Hamiltonian 
describes the mixing between d and / electrons. It was recognized for a long time 
that the d-f mixing is very im portant to u nderstand electronic and magnetic prop- 
erties of the Eu chalcogenides! I I I I So we have taken into account the d-f 
mixing term in the Hamiltonian. 

The Hamiltonian is diagonalized by the Bogoliubov transformation and the 
Bogoliubov- de Genne s equation is numerically solved by the recursive Green func- 
tion method! I I I We calculate the Matsubara Green function, 

where iOn = (2n + \)tiT is the Matsubara frequency, n is an integer number, and T 
is a temperature. The Josephson current is given by 

w 
Ij{cf>) = -tetTj2 E Tr [G.Ar',r) ~ G^Ar,r')] (9) 

UJn 7n—l 

with r' = r + a;. In this paper, 2x2 and 4x4 matrices are indicated by • " • and • ~ •, 
respectively. Throughout this paper we fix the following parameters: W = 25, /i = 
2i, and Aq — O.Oli, T — O.OlTc {Tc is the superconductor transition temperature). 



4. Josephson current through the fully polarized ferromagnetic 
insulators 

We first discuss the Josephson current through a FPFI [Fig. l(a)].'^ The phase 
diagram depending on the strength of 14x ( < Vcyi/t < 8 for FM and Vcx/t > 8 for 
FI) and Lp is shown in Fig. [3] The black (white) regime corresponds to the tt (0) 
junction [Ij — —{+)Icsm<j)]. In the case of FPFI, the tt junction can be formed. 
Moreover, the O-tt transition is induced by increasing the thickness of the FI barrier 
Lp- More detailed discussion and a physical origin of above results will be given in 
elsewhere. -^ 



5. Josephson current through the Eu-chalcogenides 

In this sect iori, we consider the Josephson transport through the Eu- 
chalcog enides.l^^^^In ca lculation, we use the following parameters in consideration 
of EuOpIlSHIMiQllTIia i^ ^ 1.25eV, g = gd + .g/ = 1.12eV, t/ = 0.125eV, and 
V^i = 0.528eV. 

We first discuss the Josephson current through the spin-filtering barrier only, 
i.e., the d-band [Fig. 4(a)]. The phase diagram depending on the strength of V^ ( 
< V^^/td < 6) and the thickness of FI (Lp) is plotted in Fig. 4(b). In this case, the 
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Fig. 3. The phase diagram depending on the strength of Vex and Lp for the FM (0 < Vcx/t < 8) 
and the fully polarized FI (Vix/i > 8). The black and white regime correspond to the it and 
junction, respectively. 
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Fig. 4. (a) The density of states for each spin direction for the spin-filtering barrier {5d band of 
Eu). (b) The phase diagram depending on the strength of Vex and Lp for the the spin- filtering 
barrier. In this case, no n junction is formed. 



TT junction is not formed irrespective of Lp and Vc^ 
dose not lead to the 7r-junction behaviors. 



Therefore, the spin filter-effect 
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Next we consider the Josephson transport through the Eu-chalcogenides with 
both d and /-bands. In calculation we set Lp — 5 and systematically change the 
values of the exchange splitting of / bands V/^{= 0.0 ^ lO.OeV) and the d-f hy- 
bridization Vdf{= —1.25 ^ 1.25 eV). Fig. 5 shows the phase diagram that was 
numerically obtained. The n junction can be realized at the certain values of Vdf 
and V/^. We found that the tt junction is formed if (1) d and / bands are over- 
lapped each other and (2) the d-f hybridization Vdf is strong enough. More detailed 
discussion for above results will be given in elsewhere."^^ 

6. Summary 

To summarize, we have studied the Josephson effect in S-FI-S junction by use 
of the recursive Green's function method. We found that n junction and the O-tt 
transition is realized in the case of FPFI. On the other hand, in the case of the Eu 
chalcogenides, the tt junction can be formed if the d and / bands are overlapped and 
the d-f hybridization is strong. Such FI based tt junctions may becomes a element 
in the architecture of " quiet qubit" . 
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